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Abstract— In this paper, an oil-heating tank is 
presented, developing its respective modeling through 
the application of the laws of conservation of mass and 
energy. Then a set of nonlinear differential equations 
are obtained representative of the process. 
Subsequently, this model is linearized obtaining its 
linearized model. With this linear model, the design of 
two control systems is presented, based on a PI 
feedback output compensator and based on a linear- 
quadratic regulator, where the control objectives are 
the storage temperature and the fuel oil level. Finally, 
the process is simulated, verifying the feasibility of the 
operation with the designed control system. 
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I. INTRODUCTION 


In most of the industries that handle fuel oil, it is 
required to store this product in tanks for later use. 
Precisely the storage condition, can generate that the 
physical properties of the fuel oil can change 
undesirably, as its viscosity, for example, due to the 
temperature gradients that could be affected. 
Therefore, a temperature regulation system is 
required and necessary to maintain its properties 
undisturbed. 


An oil-heating tank is the simple and safe 
solution for maintaining the temperature and thus 
preserving the viscosity of the fuel oil. Therefore, in 
addition to having the necessary tank for storage, it 
is also required to have a thermal system to supply 
energy to the fuel oil. It is therefore required to have 
a temperature control system. In addition, as the tank 
has overflow safety levels, a level control is also a 
must. 


It is interesting to note that there are not many 
works that address these two control goals, in fact, 
in most of them, they cover the temperature control 
as the main one, leaving in second place, the level of 
the tank, assuming steady state conditions regarding 
the level [1]-[5]. 


In this paper, an oil-heating tank is presented, 
developing its respective modeling through the 
application of the laws of conservation of mass and 
energy. Then a set of nonlinear differential equations 


are obtained representative of the process. 
Subsequently, this model is linearized obtaining its 
linearized model. With this linear model, the design 
of two control systems is presented, based on a PI 
feedback output compensator and based on a linear- 
quadratic regulator, where the control targets are the 
storage temperature and the fuel oil level. Finally, 
the process is simulated, verifying the feasibility of 
the operation with the designed control system. 


II. PROCESS DESCRIPTION 


The process and instrumentation diagram of the 
proposed oil-heating tank is presented in Fig. 1. 
From Fig. 1 it can be seen that the system is 
configured with a tank named TK-100 which is 
working as a heat exchanger. It has three inlet pipes 
and one outlet pipe. The first inlet is the volumetric 
flow of the oil to be processed, which flows from an 
external process unit. This inlet is associated with 
the following dynamic variables fi(), 7i(f), and 
P,(0), which are the flow, the temperature, and the 
pressure, respectively. In addition, the piping of this 
inlet is connected to the tank at a height of Hy 
measured form the ground. The second inlet is the 
saturated steam flow, supplied by a generator unit 
external to this process, which allows to supply heat 
to the oil, through its flow, by a heating coil located 
inside the tank. The heating coil is built by a metal 
tube of thickness 20 BWG. The dynamic variables 
of the steam are the mass flow w,(A), the temperature 
T,(¢), and its saturation pressure P,(¢). The third input 
is the constant volumetric flow of inert gas fy which 
allows maintaining a constant pressure (P2) inside 
the tank above the oil level, due to the presence of a 
blanket of this inert gas. Finally, the output of the 
TK-100 is the oil already processed and reaching the 
reference temperature (7°3(4)). The dynamic 
variables of the outlet are the flow /(d), its 
temperature 73(t), and its pressure P3(?). 


The process assumptions are that the TK-100 is 
well insulated, the physical properties of the oil are 
not relatively strong temperature-dependent 
functions [6], [7], the liquid (oil) inside the tank is 
well mixed, the liquid level completely covers the 
heating coil, constant and uniform densities, specific 
heats, and steam latent heat. Moreover, the metal of 
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Fig. 1. Oil-heating tank process. 


the tank and heating coil is at the same temperature 
as the condensing steam, negligible heat losses and 
transportation lags, and finally, constant pressure 
drop across the steam valve. Condensation from the 
heating coil is trapped by a vapor trap (trap) and 
conducted to another process unit, also external. 


From the above assumptions, it can be concluded 
that the temperature of the oil inside the tank is the 
same as that of the outlet, i.e. 73(4). The height of the 
level related to the liquid inside the TK-100 is 
labeled as h(?). 


The volumetric flows, mass flow, temperatures, 
height, and pressures are measured in m?/s, kg/s, °C, 
m, and kPa respectively. 


From Fig. 1 it can be inferred that the control 
philosophy corresponds to controlling 73(¢) and h(d), 
which are the controlled variables of the system, 
through the manipulation of valves TV-101 and LV- 
102 respectively. The temperature 73(f) and the 
height A(f) are measured through the transmitters 
TT-101 and LT-102 respectively. Finally, the 
controllers associated to the controlled variables are 
TIC-101 and LIC-102 respectively. On the other 
hand, there are two uncontrolled valves whose 
position is fixed. These valves are HV-104 and HV- 
103 and throttle the flows fy and f3(4) respectively. 

Ill. PROCESS MODELING 

The mass and energy balances [6]-[8] developed 
for the process shown in Fig. | are presented below. 
From the synthesis of these balances, the nonlinear 
model of the system is obtained and shown. 
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A. Mass and Energy Balance 


Through the use of the law of conservation of 
mass, the differential equation that relates the height 
A(f) is obtained and presented in (1). 
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where ATK-100 is the cross-sectional areas of the 
TK-100 in m2. This area is uniform throughout. 


On the other hand, the energy balances around 
the TK-100, which is related to 73(4) and the heating 
coil, which is also related to 7,(#) are presented in (2) 
and (3) respectively. It should be noted that for 
deriving the energy balance related to the TK-100, it 
is assumed that the volume of the liquid in the tank 
is fairly constant, for a given operating height, i.e., 
the level control LIC-102 operates properly. 
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where p, cy, and cy are the density and heat capacities 
at constant volume and pressure of the oil, measured 
in kg/m? and in kJ/kg-°C respectively. In addition, 
U, Vrx-100, An, Cm, and As are the total heat transfer 
coefficient, the effective volume (= Arx-100-/) of the 
tank, the cross section and heat capacity of the 
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heating coil, and the latent heat of condensation 
respectively. Each of these process constants are 
measured in W/m?-°C, in m%, in m’, in kJ/°C, and in 
kJ/kg respectively. 


From Fig. 1 and according to [1] the expressions 
of the flows w,(t), fi(A, and f3(t) are a function of the 
pressure drops, through the valves TV-101, LV-102, 
and HV-103 respectively and are presented in 
equations (4)-(6) (it is not necessary to consider the 
dynamics of the flow fy, since this flow does not play 
a role in the balances in (1)—(3). 


w,(t) = Cry-101 * VPrv-101(t) (4) 
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Here g and G; are the gravitational acceleration and 
specific density of the oil measured in m/s? and 
dimensionless, respectively. Moreover, Crv-101, CLv- 
102, and Cyy-103 are the valve coefficients of the TV- 
101, LV-102, and HV-103, respectively. 


Replacing (4)-(6) in (1)-(3) the non-linear 
model of the oil-heating tank is obtained and 
presented in (7). 


From (7) it should be noted that Cu = Pm:Lm:Cpm 
where Om, Lm, and Cpm are the mass, linear length, 
and heat transfer coefficient of the metal tube of the 
heating coil [6]. The units of these constants are 
kg/m, m, and kJ/kg-°C respectively. 


In the following section the linear model of the 
system will be obtained and presented. 


B.  State-Space Model of the Process 


In order to design a linear controller, such as a 
proportional-integral (PI) compensator or a linear- 
quadratic regulator (LQR) [9]-[11], it is required to 
linearize the model of (7), thus obtaining the state- 
space model. The linearization obtained is based on 
the steady-state operating points [12], [13] that are 


(As-crv—101°VPTv-101 (t)-U-Ap-(Ts(t)-T3 (t))). (7) 


the outcome of developing the system of equations 
shown in (8) and that is achieved by forcing the 
derivatives to zero, ie., dh(t)/dt = d73(H/dt = 
d7,(0)/dt = 0. Then, replacing the time variables by 
their steady-state variables. Note that the steady- 
state variables are presented by capital letters and a 
superscript ss. 


The unknowns of the system of equations in (8) 
are the positions of the valves LV-102 and TV-101, 
Le., VPS~y-102 and VP*ty-101 respectively, and the 
saturated steam temperature, ie., T*;. The other 
variables are all known and given. Table I lists all 
the constants and steady state variables of the 
process. The solution of (8) are presented in 
equations (9)—(1 1). 


Having the steady state operating points it is 
possible to linearize the model in (7) [12]. Applying 
the Taylor series expansion to (7) around the 
operating points in (9)—(11), the linear model, as a 
function of deviation variables, of the process is 
obtained and shown in (12). In this case, where all 
auxiliary constants present in (12) are defined in 


(13)—(15). 


The deviation variables are A(t) = h(t) — H*S, 
T3(t) = T3(t) —T3"*, 7,@) =T.@) — TS, 
OPv-102(t) = vPiv-102(t) — VPiv-102, Pi) = 
Pi(Q)- PP, PO =P, 10) =) - Tr, 
PPrv-101(t) = vPry-101(€) — VPpy_101. Finally, 
from the linear model of the process described in 
(12), the state-space model is obtained and described 
in (16). 


The vectors and matrices are defined in (17) and 
(18) respectively [12]-{14]. Here {%, 7} < {R?} and 
fie {R°}, {Am, Cm} € x3 {R}, and {Bm, Dm} € 
AMxs {R}. 


IV. CONTROL SYSTEM DESIGN 


The following section presents the designs of 
two linear controllers, which are applied to the oil- 
heating tank. The first controller is based on a PI 
compensator and the second one is based on an 
LQR. 
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TABLE I. CONSTANTS AND STEADY STATE VARIABLES OF THE 


PROCESS 
Variables Values 
p 848.979 kg/m? 
Cp 1.8841 kJ/kg-°C 
Cy 1.8841 kJ/kg-°C 
Gr /1000 = 0.8490 
g 9.81 m/s? 
An 11.845 mm? 
La 296.875 m 
Pm 0.265 kg/m* 
Com 0.502 kJ/kg-°C 
U 771 W/m?-°C 
Cu 39.4933 kJ/°C 
As 2,030.6 kJ/kg 
D (tank diameter) 3m 
A 0.25-1-D? = 7.0686 
TK-100 2 
m 
max (Maximum tank height) 2.5m 
Vrx-100 7.0686 m? 
Ay 1.5m 
HY 0.5-Amax ¥ 1.25 m 
Ps, 10 m/s 
T*; 21.11 °C 
T°; 93.33 °C 
Te: 175°C 
PS. 792.897 kPa 
PP; 31.026 kPa 
PS, 276 kPa 
PS, 103.421 kPa 
VP Lv-102 0.31 per unit 
VP*ty-101 0.041 per unit 
CLV-102 8.57-10° 
CHV-103 8-10° 
CTV-101 1.174 


The data of the transmitter are the gains and time 
constants of each transmitter, 1.e., Krr-io1, Kir-102, 
Trr-101, and 7-102 respectively are listed in Table II. 
In addition, the dynamics of the transmitters are 


modeled as first order responses [1]. Although the 
TT-101 does not perform any conversion in the 
measurement, i.e., it reads degrees °C, the LT-102 
converts the measurement from m to %, since, in 
normal, tank capacities are measured as a percentage 
of their maximum height, i.e., hmax. The scale of 
TT-101 and LT-102 are [0, 3.048] m and [0, 149] °C 
respectively. 


A. Pl-Compensator Design 


By applying the Laplace transform to the linear 
model (12), considering zero initial conditions, the 
open-loop model of the oil-heating tank is shown in 
(19). In addition, the block diagram of the linearized 
process is depicted in Fig. 2. 


From (19), s is the Laplace operator. From Fig. 2 
it is possible to identify the two transfer functions 
related to the control loops TIC-101 and LIC-102. 
Regarding TIC-101 the transfer function relating 
T3(s) and VPry_101(s) is presented in (20). Related 
to LIC-102 the transfer function relating H(s) and 
VP.y-102(s) is shown in (21). 


A(s) = K12 
VPLy-1026) = s—Ky1 (20) 
T3(s) = K21'K31 (21) 


VPpy—101(9) — (s—Ka1)-(s—K31)" 


In order to design the PI compensators for TIC- 
101 and LIC 102, the transient analysis (step 
response) and root locus diagrams of (20) and (21) 
are obtained and shown in Fig. 3(a) and (b), and Fig. 
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4(a) and (b) respectively. From Fig. 3(a) and 4(a) it 
can be seen that both A(t) and T3(t) reach a steady 
state value, where the slowest dynamics is that 
corresponding to 7;(t). On the other hand, and 
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according to Fig. 3(b) and Fig. 4(b) it can be 
concluded that the transfer functions shown in (20) 
and (21) are stable, since all the poles are located on 
the left side of the s-plane. 
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TABLE II. PARAMETERS OF TRANSMITTERS TT-101 AND LT-102 


Parameters Values 
Krr-i01 1 TO%/°C 
Kwr-102 109.769 TO%/m 
Trr-101 0.05 s 
1-102 0.01 s 


TABLE III. TUNING PARAMETERS OF COMPENSATOR Gei(S) AND 


G.2(S). 
Tuning parameters Values 
Koct 0.006 
Kit 9.152-10° 
Kyc2 0.096 
Kier 3.96-104 


Using this information provided by Fig. 3 and 4 
and with the aid of MATLAB-Simulink (with the 
SISOTOOL command), the PI compensators related 
to the TIC-101 and LIC-102 control loops are 
designed. It should be noted that SISOTOOL is 
based on the root locus analysis method [12]-[14]. 


Fig. 5 shows the block diagram of the closed- 
loop process, including the PI compensators and the 
TT-101 and LT-102 transmitters (red-dash lines). 
Fig. 5 shows the saturators connected to the output 
of the PI compensators, in order to prevent the 
control valves, i.e., TV-101 and LV-102, to receive 
non-binary signals. Finally, the transfer function of 
the compensators is shown in (22). Here kp and kj are 
the proportional and integral gains of the G 
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The design of the PI compensators related to the 
LIC-102 and TIC-101 controllers, i.e., Goi(s) and 
G.2(s) respectively, provides the tuning parameters 
of each of the compensators listed in Table IT. Fig. 
6(a) and (b) show the responses to a step change 
associated with the TIC101 and LIC-102 controllers, 
respectively, which are described by the red dash 
lines in Fig. 5. From Fig. 6, it can be observed that 
both measured variables, i.e., T;(t) and A(t) reach a 
steady state value, showing that the dynamics of 
LIC-102 and TIC-101 are stable. 


B. LOR Design 


Regarding the design of the LQR, which is based 
on optimal control techniques [9], [10], the state 
equation of (16) is considered, i.e, = Ay *X®+ 
By * Hi. On the other hand, the optimal control law 
ii = —K -X is defined in order to minimize the cost 
function shown in (23). 
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where Q and R are positive definite Hermitian 
matrices [9], [10]. Matrix K is a matrix with optimal 
values that is determined from the solution of the 
Riccati equation [9], [10], [12] and matrices Q and 
R are defined in (24). 


As the LQR algorithm is conceptualized as a 
regulator and thus operates with zero references [9], 
[10]. it is necessary to customize the algorithm to 
apply it to the oil-heating tank, since, in this case, 
there are references different from zero, obviously. 
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Fig. 4. Transient analysis and root locus diagram regarding (21). (a) Step response. (b) Root locus diagram. 
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Fig. 5. Closed-loop block diagram of the oil-heating tank model, including PI compensators and measurement transmitters. 
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Fig. 6. Step response of closed-loop transfer functions LIC-102 and TIC-101. (a) TIC-101. (b) LIC-102. 
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First, the steady-state vector is defined, which 
operates as a vector that takes into account the 
reference values, i.e. X° = [H®, 73, T,S]7, X € 
{R3}. A vector is also defined based on the steady 
state inputs, ie. USS = [VPtyv-102, VPrv-i01, Pris, 
P38, T)5]', US € {R>}. Then, redefining the dynamic 
equation in (16) and the control law, in terms of the 
absolute and steady state variables (see subsection 
3), the optimized model involving the vectors X* 
and U® is obtained and shown in (25). In this case x 
=[h(), T3(), TAO)", x € {R°}, w= [vpiv-102(), vprv- 
oi(Z), Pi(t), P(t), T(]", w € {R°}. 


By using the MATLAB command jgr, the 
optimal matrix K is obtained and shown in equation 
(26). The eigenvalues of the optimized model (25) 
are presented in vector E (26), which correspond to 
negative values, indicating that (25) is stable [9], 
[10]. Fig. 7 shows the LQR diagram incorporated in 
model (25). 


(oo (25) 
(u — U*S) = —K: (x— X*$) 

4.8344 2.7554 0.0038 

—8.8926 1.5386 0.0035 
K=]| 0.0000 0.0000 0.0000], 

—0.0002 0.0000 0.0000 

—0.0533 0.0092 0.0000 

—0.0924 
—231.6433 


V. SIMULATION RESULTS 


Simulations of the lossless linear model 
represented by (16) were carried out in MATLAB- 
Simulink. The process parameters and variables are 
listed in Table I. 


It is assumed that the oil-heating tank is 
operating at steady state conditions with a level of 
A = 50 % with an oil temperature inside the TK- 
100, 73 = 93.33 °C, a saturated steam temperature, 
T*; = 175 °C and a pressure P*; = 31.026 kPa. Then 
at 800 s there is a step change in h*(f) reaching a new 
steady state value of 30 %. Then at 1,000 s, Pi(f) 
under-goes a perturbation decreasing the inlet oil 
pressure to 26.3721 kPa, and finally at 1,200 s a step 
change in 73(¢) is generated reaching a new steady 
state value of 70 °C. 


Fig. 8 shows the dynamics of the variables h(d), 
T3(0), vptv-i0i(0), vptv-102(2), and P\(0), when the 
controller corresponds to a PI compensator and a 
LQR. Fig. 8(a) shows a comparison of the dynamics 
of h(f), when the controller is an LQR (the blue one) 
and when it is a PI compensator (the cyan one), in 
addition the dynamics of h(t) before being measured 
by the LT-102 transmitter is displayed (the green 
one). The red one is the reference value. Fig. 8(b) 
shows the comparison of the dynamics of 73(¢) in the 
same format as Fig. 8(a). Fig. 8(c) shows the 
comparison of P;(¢) dynamics and finally Fig. 8(d) 
and (e) are the comparisons of the dynamics of vptv- 
102(¢) and vprv-101(f) valve positions respectively. It 
can be seen from Fig. 8(a) and (b) that the shortest 
settling times are presented by the dynamics /(f) and 
T3(t) when the controller is an LQR instead of a PI 
compensator. In addition, from these same figures it 
is clear to see the effect that the time constants of the 
transmitters have on the display of the dynamics. 
The green dynamics shows the initial state of the 
level H** (50 %) and temperature 73° (93.33 °C), 
instead of the cyan figure, which is the dynamics of 
the level and temperature, coming from the 
transmitter LT-102 and TT-101 respectively, which, 
due to the effect of the time constants, the initial state 
of h() and 73() are shown as 0 %. 


On the other hand, from Fig. 8(a) and (b) it can 
be seen that when TIC-101 and LIC-102 controllers 
are PI compensators the step change in h'(t) and 
Ts (0) generates disturbances in less magnitude than 
when the controllers are combined in a LQR in the 
variables 73(¢) and h(t) respectively. However, as 
mentioned above, the settling time is much longer 
when PI compensators are involved. Moreover, Fig. 
8(c) also shows that the temperature dynamics 7,(4), 
which is measured at the heating-coil output (see 
Fig. 1), has more stable dynamics when an LQR is 
applied instead of PI compensators. This feature is 
of key importance, since it is essential in this type of 
process to keep the temperature of the saturated 
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Fig. 8. Simulation results under transient operation. Reference level h’(¢) = 50 %, temperature 7” ;(f) = 93.33 C, and steady-state value of 
P,(t) = 31.026 kPa. A step change in h’(¢) at 800 s, a step-change in 7°;(¢) at 1,200 s, and a step-change in P,(¢) at 1,000 s. (a) Level 
dynamics h(t). (b) Temperature dynamics 7;(t). (c) Temperature dynamics T,(¢). (d) Pressure dynamics P(t). (e) Valve position dynamics 
vprv-102(t). (f) Valve position dynamics vpry-ioi(¢). All dynamics are overlapped and match when the controller is a LQR and when it is a 


PI compensator. 


steam as stable as possible in order to reduce the cost 
of saturated steam production. It is well known that 
the generation of saturated steam is one of the main 
expenses in thermal or petro-chemical plants, where 
there are a large number of oil-heating tanks. It is 
interesting to note that, when an external disturbance 
occurs in P,(f) (see Fig. 8(d)), none of the variables 
presented in Fig. 8 is affected, showing then that the 
process with both types of controllers (PI and LQR) 
presents an important robustness. Finally, from Fig. 
8(e) and (f) it can be seen that the behaviors of the 
valve positions vpty-ioi(4) and vpry-102(t) present a 
smoother modulation when the controllers are PI 
than when it is a LOR, However, when the controller 
is an LOR, vprv-101(¢) and vprv-102(f) are more stable 


but involve some peaks, resulting in openings and 
closings at 100 % and 0 % respectively in a very 
short time, which could be harmful to the 
mechanical elements of these valves. 
VI. CONCLUSION 

The study of an oil-heating tank incorporating 
mathematical modeling and the de-sign of two linear 
controllers has been presented in this work. Based 
on the law of conservation of mass and energy 
applied to this process, the nonlinear dynamic 
equations that model the oil-heating tank are 
obtained. Then, the operating points of the steady- 
state process are calculated to finally linearize the 


system using the Taylor expansion around the 
calculated operating points. As a result of the 
linearization, a state-space model of the oil-heating 
tank is obtained. 


Using the state-space model, two proportional- 
integral compensators are designed; one to control 
the level of the oil-heating tank /(A) and the other to 
regulate the outlet temperature of the processed oil 
T;(t). On the other hand, to regulate these same 
variables, a linear quadratic regulator is designed 
and the results of these controllers are compared. 


In overall, and according to Fig. 8, it can be seen 
that the process with both types of linear controllers, 
i.e., the proportional-integral compensators and the 
linear quadratic regulator works properly, but in 
particular it can be concluded that the dynamics of 
the controlled variables, i.e., A() and 73(¢) when the 
process operates with two proportion-al-integral 
compensators working as controllers, present a 
longer settling time compared to the case when the 
process operates with a linear quadratic regulator. In 
addition, the dynamics of h() and 73(4) show a better 
smoothness against produced disturbances. 


On the other hand, when the process operates 
with a linear quadratic controller, the dynamics of 
Ts(¢) is much smoother than when operating with 
proportional-integral compensators, being a key 
factor the stability of this temperature, since the cost 
of the production of saturated steam has as a 
function, precisely the stability of its temperature. 


Finally, it can be concluded that when the 
process is controlled with the proportion-al-integral 
compensators, the dynamics presented by the valve 
positions vprtv-ioi(f) and vpty-1o2(t) present a 
smoother modulation although more variable than 
those shown when the controller is a linear quadratic 
regulator, However, when the process is controlled 
by the latter, the valve positions show certain peaks 
of 100% opening and absolute closures in very short 
periods of time, which can cause important damages 
to the mechanical elements of the valves. 
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